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ABSTRACT 

We investigate whether the Hubble sequence can be reproduced by the relics of merger events. We verify that, at Zmedian=0 .65, the 
abundant population of anomalous starbursts - i.e. with peculiar morphologies and abnormal kinematics - is mainly linked to the local 
spirals. Their morphologies are dominated by young stars and are intimately related to their ionised-gas kinematics. We show that 
both morphologies and kinematics can be reproduced by using gas modelling from Barnes' (2002) study of major mergers. Their gas 
content may be indirectly evaluated by assuming that distant starbursts follow the Kennicutt-Schmidt relation: the median gas fraction 
is found to be 31%. Using our modelling to estimate the gas-to-stars transformation during a merger, we identify the gas fraction in 
the progenitors to be generally above 50%. 

All distant and massive starbursts can be distributed along a temporal sequence from the first passage to the nuclei fusion and then to 
the disk rebuilding phase. This later phase has been recently illustrated for J033245. 1 1-274724.0, a distant compact galaxy dominated 
by a red, dust-enshrouded disk. This active production of rebuilt disks is in excellent agreement with model predictions for gaseous 
rich encounters. It confirms that the rebuilding spiral disk scenario - a strong and recent reprocessing of most disks by major mergers 

- is possibly an important channel for the formation of present-day disks in grand-design spirals. Because half of the present-day 
spirals had peculiar morphologies and anomalous kinematics at Zmedimi=0-(>5, they could indeed be in major mergers phases 6 Gyrs 
ago, and almost all at z~ 1 . 

It is time now to study in detail the formation of spiral disks and of their substructures, including bulge, disks, arms, bars, rings that 
may mainly originate from instabilities created during the last major merger. Many galaxies also show an helicoidal structure, which 
is probably due to a central torque, and seems to play an important role in regulating the angular momentum of the newly-formed 
disks. 

Key words. Galaxies: formation -Galaxies: spiral - Galaxies: kinematics and dynamics 
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1. Introduction 



The tidal torque theory (Peebles ' 1 976 '; White' 1 984") assumes that 
the angular momentum of disk galaxies had been acquired by 
early interactions. This theory has been supported for a long 
time: in fact almost all massive galaxies are regular, including 
rotational disks and dispersion-supported bulges or their mixes, 
and they outline the local Hubble sequence. However if spiral 
disks were formed at early epochs - z > 2 - they could have en- 
countered severe damages from later major interactions. Galaxy 
collisions appear to be too frequent, to allow many disks to sur- 
vive (.Toth & Ostriker 1992j , and this might happen even at z<l 
( |Hanimer et al.|2009a| . 

How was the Hubble sequence 6 Gyrs ago? Galaxy mor- 
phologies s trongly evolve (|van den Bergh 2002 2009 Zheng 
|et al.|2004{ |Conselice et al.|2005| ) towards much more peculiar 
structures. This combined with the coeval evolution of star for- 
mation rate and stellar mass densities, of O/H gas abundances 
and pair statistics, has prompted us to propose the disk rebuild- 
ing scenario (Hammer et al. 2005j). This scenario describes most 
of these evolutions as due to a recent merger origin of most spi- 
rals. To reproduce the observed evolution requires that 50 to 75% 
of the present-day spirals have been formed - i.e. their disks re- 
processed by mergers - during the last 8 Gyrs (z< 1). Within 
this period, observations tell us that more than half of the stel- 



lar mass in spirals has been formed and this can happen through 
gas compression occurring during the different phases of major, 
gas-rich mergers (Hammer et al. 2005 1. In this theory, the disk 
angular momentum mostly results from the orbital angular mo- 
mentum of the last major collision ( jPuech et al.|2007at|Hammer| 
et al. '2007 ; Hopkins et al. 2009a). A more recent epoch for disk 
formation is indeed supported by the large decrease with red- 
shift of the fraction of rotationally supported disks. |Neichel et aL] 
( |2008| l (hereafter IMAGES-II) found that rotational disks were 
two times less abundant at Zmedian=0.65, a result that is based on 
a study combining detailed-morphology and spatially-resolved 
kinematics. 

How can mergers be related to the regular local galaxies, 
of which our Milky Way has been so often taken as typical? 
Deep observations attest the rather tumultuous history of sev- 
eral nearby galaxies that is imprinted in their inner halo (e.g. 
M31, see [Brown et al.|2008|pata et al.|20()5| see also |Davidge" 



2008] for M81). The Milky Way appears to be quite exceptional 



(Hammer et al. 2007), possibly related to its quiescent merger 
history. Particularly, its halo is the most primordial within galax- 
ies with similar masses ( [Mouhcirie 2006 1 and it shows an angular 
momentum two times smaller than that of typical spirals. 

Galaxy simulations can help to test various galaxy forma- 
tion scenarios. Assuming large accretions of cold gas flows may 
reproduce several correlations, mostly those linking the gas con- 
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sumption and the assembly of the stellar mass ( [Dekel et al.| 
2009| l. There is, however, no convincing observation of signif- 
icant cold gas-flow in local or distant galaxies, while mergers 
are well identified in the local and distant Universe. Ho pkins et| 
al. ( 2009a[ l (see also Robertson et al.|[2006 ) successfully tested 
disk survival during merging. Resulting mergers are producing 
disks whose angular momentum differs in direction and ampli- 
tude from those of the progenitors. The predictions of H opkins, 
|et al . (2009a) rely on processes dominated by pure gravitation, 
which are indirectly affected by feedback effects. This is cer- 
tainly true for massive galaxies for which stellar feedback - 
including outflows- is unlikely to be an efficient way to redis- 
tribute the material during a collision. If the gas fraction is suf- 
ficient (about 50%), they predict that the re-formed disk can be 
the dominant component in the reshaped galaxy. 

The rebuilding disk scenario proposes a merger origin for 
spirals and by extension of the whole Hubble sequence, from 
ellipticals (Toomre &Toomre ||1972[ |Toomre||198l] [Barnes &| 



Table 1. Morpho-kinematical classification of 52 

dSOOSl) (see their 



Table 4); for 



Hernquist 



1992| i to late type spirals. In fact, the orbital angu- 



lar momentum provided by major mergers could solve the an- 
gular momentum problem (Mailer, Dekel, '& Somerville|2002| l. 
Considerable work is, however, needed to support the scenario. 
Could a cosmological distribution of orbital parameters rebuild 
small bulges within some rebuilt disks, and large bulges within 
others? Could it reproduce the Hubble sequence statistics of 
bulge to disk ratio? Hopkins et al. ( [2009a b ) partially brought 
a positive answer to this question. Their model indeed recov- 
ers the well-known correlation between bulge-to-disk ratio and 
mass. 

Observations may prove or invalidate the rebuilding disk sce- 
nario. It is well known that the gaseous content of galaxies in- 
creases rapidly with redshift. But does it reach the values re- 
quired to rebuild a disk in case of mergers? We must also ob- 
serve the details of the physical processes in galaxies at differ- 
ent epochs and examine directly their evolution. At very high 
redshift (z> 2) cosmological dimming prohibits the examination 
of the optical radius of a disk even with the largest space tele- 
scopes, and it is difficult yet to gather a representative sample 
of such galaxies. At intermediate redshifts the situation is much 
better even if it needs pushing the present observational tools 
near their limits. Up to z-QA and z=1.3 the optical radius of a 
redshifted Milky Way can be retrieved from GOODS and UDF 
imageries, respectively. 

The IMAGES study aims to identify the physical processes 
that link distant (z ~ 0.65) to local galaxies. Its selection is lim- 
ited by absolute J-band magnitude (My(AB)< -20.3), a quantity 

hereafter 



2008 



relatively well linked to the mass (Yang et al. 
IMAGES-I). Using such a limit, [Delgado-Serrano eTall ( [2009] l 
have shown that z ~ 0.65 galaxies have to be the progenitors 
of local galaxies selected in a similar way. They find that the 
fraction of E/SO did not evolve since the last 6 Gyrs, while spi- 
ral galaxies were 2.3 times less abundant. They indeed use a 
quite restrictive method to classify morphologies, assuming that 
spiral galaxies in the past should have similar properties than 
what they possess today. Using such a morphological classifica- 
tion, IMAGES-II ( |Neichel et"aL][2008 ) demonstrated an excel- 
lent agreement between morphological and kinematical classifi- 
cations. In other words, most rotating galaxies (80%) show spiral 
morphologies while most galaxies (90%) with anomalous kine- 
matics present peculiar morphologies. The above results have 
an important impact: anomalous kinematics of the gaseous com- 
ponent (from the ionised gas, [(9//]/13726;3729 lines) is almost 
always linked to anomalous morphological distribution of the 
stars. Altogether the above results imply that more than half 



galaxies from Neichel et al. 

comparison, the last column shows the fractions derived from 

for galaxies in the same 



the SDSS ( [Nakamura e"raL][2004 



mass range (e.g. Hammer et al., 2005). CFRS (Hammer et al 



11997 1 found that 60% of Zmedian=0-65 galaxies have spectra with 
and are classified as starbursts. References are 



Woi OII) >15A 

■ Neichel et al.| ( [2008) : N08, .Zheng et al.,(2004) : Z06, ,Nakamura| 
|etal.[ ( |2004| i:N04 and [Hanuner et al.[ ( |2005| l: H05. 
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of the present-day spirals had anomalous kinematics and mor- 
phologies, 6 Gyrs ago. 

Anomalous galaxies are also responsible ( Flore s et al.|2006| l 
of the most striking evolution of the TuUy-Fisher relation, i.e. 
it is heavily scattered at z~0.6-l ( Conselice et al.|20 05). Major 
mergers can reproduce this evolution ( [Covington et al. 2009 ) as 



well as a similar trend for the jdisk-yjiai relationship (Puech et al. 



|2007a I. The goal of this paper is to verify whether the observed 
evolution can be mostly related to merger events, i.e. to test if the 
rebuilding disk scenario is consistent with the observed evolu- 
tion of morphology and kinematics. We thus defined 3 different 
morpho-kinematical classes following the Table 4 of [Neichel et[ 
[ai:i( |2008] l: 

- rotating spiral disks are galaxies possessing a rotating velocity 
field, including a dispersion peak at the dynamical centre (see 
e.g. Flores et al.|2006) , and showing the appearance of a spiral 
galaxy; 

- non-relaxed systems are galaxies that have a strong discrep- 
ancy from a rotational velocity field and whose morphology is 
peculiar; 

- and semi -relaxed systems that possess either a rotational veloc- 
ity field and a peculiar morphology or a velocity field discrepant 
from a rotation and a spiral morphology. 

Table 1 summarises the statistics at Zmedian-0-65 and compares 
them to local galaxies from SDSS (Nakamura et al. 2004). Note 
that similar statistics combining kinematics and morphology 
does not exist for both local galaxies and quiescent distant galax- 
ies, and thus the corresponding fraction of non-relaxed systems 
is still inaccurate. Indeed Delgado-Serrano et al.[ ([2009]) found 
that 10% and 25% of local and quiescent distant galaxies show 
peculiar morphologies, respectively. Given this. Table 1 provides 
only a lower limit of the fraction of distant galaxies that show 
anomalous properties and that are the progenitors of present-day 
spirals. 

Morphological and kinematics properties are coming from 
the IMAGES survey and the complete picture has been also 
provided by deep photometric and spectroscopic measurements 
necessary to estimate their SFR (both in UV and IR, see also 
Puech et al.|2009a ), their stellar masses and their chemical and 
stellar population decomposition. Section 2 describes our pro- 
cedure to test whether or not morpho-kinematical properties of 
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distant starbursts can be reproduced by galaxy interactions pro- 
cesses, including during the remnant phase. In section 3 we show 
the overall properties of distant starbursts, including their gas 
richness that is of crucial relevance to infer whether mergers 
may lead to disk rebuilding. In section 4, we discuss the results 
and conclude in section 5 on the validity of the disk rebuilding 
scenario. Throughout the paper, we adopt Ho - 70 km/s/Mpc, 
Q.M — 0.3, Q.A — 0.7 and the AB magnitude system. 



2. Could distant starbursts properties be 
reproduced by merger or merger remnants? 

2.1. Detailed analyses of individual distant galaxies 

Detailed analyses of four distant galaxies of the IMAGES study 
have been performed by Puech et al., ( 2007a 2009a), Hammer] 
[etaLl ( |2009bl l and |Peirani et al.lpOOsj , and four other studies of 
individual galaxies are in progress (Yang et al., 2009; Fuentes 
et al. and Peirani et al. in preparation). By modelling gas mo- 
tions as well as morphologies, these studies have shown their 
ability in reproducing the properties of distant galaxies with a 



similar accuracy to what is done for nearby galaxies. Puech et 



al. (2007a) have demonstrated that spatially resolved kinemat- 



ics is sufficiently sensitive to detect the infall of a 1:18 satellite 
in a z=0.667 galaxy. Peirani et al. (20081 identified a giant and 
starbursting bar induced by a 3:1 merger, and simulated both 
morphologies and the off-centre dynamical axis. In this case, 
the gas pressured in the tidally formed bar has condensed into 
young and blue stars. Hammer et al. ( 2009b[ ) identified a com- 
pact LIRG dominated by a dust-enshrouded compact disk that 
surrounds a blue, centred helix (so-called a "two arms-plus-bar" 
structure). They interpret (see their Fig. 7) this structure as regu- 
lating the exchanges of the angular momentum and possibly sta- 
bilising the new disk ( [Hopkins et al.|2009a) l. Indeed gas inflows 
along an helix are usual in simulations of mergers, especially in 
inclined and polar orbits. This gaseous-rich galaxy appears to 
be an archetype of a disk rebuilding after a 1:1 or a 3:1 merger 
with an inclined orbit. Puech et al. (2009a) demonstrated that the 
presence of ionised gas without stars near a highly asymmetric 
disk can be only reproduced by a remnant of a merger. 

These studies have been successful because they compared 
simulations of the gas phases to observations of both the mor- 
phology and the ionised gas motions. Morphologies of star- 
bursts -especially the numerous blue or dusty regions- are mostly 
relics of gas phases recently transformed into young stars that 
ionise the gas. Thus a common physical mechanism should re- 
produce them together with the observed large-scale motions of 
the ionised gas. Within most starbursts, the light is indeed dom- 
inated by < lOOMyrs-old stars and at large distances, spatially- 
resolved kinematics only detect large-scale motions, with typical 
scales of ~3kpc. A typical motion of lOOkm/s would cross such 
a length scale during ~ 50Myrs (32 Myrs for motions parallel 
to the sky plane). Thus many morphological features with blue 
colors (bars, rings and helixes, see Peirani et al._2008, Hammer, 
et al.|2009b i should be imprints of the gas hydrodynamics and 
they can be compared to the gas kinematics. 

2.2. A general method to compare galaxy-simulations to 
distant starbursts 

For reasons of homogeneity, we study here the sub-sample 
of 33 IMAGES starbursts (see IMAGES-I) observed in the 
CDFS-GOODS. This sub-sample is representative of My(AB)< 
-20.3 starbursts (see IMAGES-I). Two galaxies have been re- 



Fig. 1. b-i-v, i and z combined images of 2 rotating disks iden- 
tified by IMAGES. Dotted line is the superposition of the dy- 
namical axis, dotted squares indicates the dispersion peak ele- 
ment. On the left, the dynamical and optical axes are aligned, 
and the dispersion peak is at the mass center, as expected for 
a rotation (see [Flores et alT][2006 l). On the right (J033238.60- 
274631.4) there is a slight misalignment of the dispersion peak 
that is likely caused by the nearby passage - 15 kpc- of a bulge- 
dominated galaxy, causing the observed burst of star formation 
at the bottom edge of the spiral galaxy. The velocity difference 
between the two galaxies is 540 km/s, a value based on spec- 
troscopy. 



jected from the original sample of IMAGES-I, one (J03 3210.76- 
274234.6) because it turns out not to be a starburst (Yang et 
al. 2009) and another one (J033250.24-274538.9) because the 
HST/ACS images are corrupted. We have verified that this sub- 
sample is representative of the stellar-mass and star formation 
densities at Zmedian-0.65 (see e.g. Ravikumar et al., |2007 1. In 
this sample we do find only 6 rotating spiral disks to which we 
add one giant spiral (J033226. 23-274222. 8 that is also rotating 
while it likely experiences a satellite infall causing a small shift 
in the observed dispersion map ( Puech et al.^2007b) . Note also 
that one of the rotating spiral galaxy (see Fig. I, right) is within 
a confirmed interaction with an elliptical galaxy. The 26 other 
galaxies all show peculiar morphologies and/or anomalous kine- 
matics and are classified as such as non or semi-relaxed systems. 

It is an Herculean task to analyse in details all the consid- 
erable amount of data for each of these galaxies, as it has been 
described for few galaxies in section 2.1. The accurate modelling 
of both morphology and kinematics takes several months, from 
two to six months for a well-experimented user This is due to 
the wide complexity of the morphologies and kinematics in these 
non-relaxed galaxies as well as the large parameter space off'ered 
by the simulations (mass ratio, orbit, temporal phase, peri-centre 
radius and parameters of the encounters, viewing angles). Our 
goal here is restricted to the following question: Can we sketch 
both the morphological and kinematics main properties of the 26 
non-(or semi-)relaxed galaxies by merger or merger remnants? 



Recently, Barnes & Hibbard (2009 ) have defined a modelling 
tool to identify merger orbital parameters. It allows to change 
many parameters including the viewing angle. However at high- 
z we cannot identify low surface brightness tidal features. We 
propose here to adapt a similar modelling tool for high-z obser- 
vations, also allowing changes of the viewing angle. We then 



used the models from Barnes (2002 ), that include 12 config- 
urations with mass ratio ranging from 1:1 to 3:1, orbits from 
INClined, DIRect, POLar and RETrograd and pericenter radii 
from rp-Q.2 to 0.4 (see ^Barnes 2002 for more details). We have 
recovered the ZENO code sourc^arid follow all the parameters 



' The ZENO simulation code was retrieved from the Josh Barnes 
website (http://ifa.hawaii.edu/~bames/software.htmll. The code has 
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in |Barnes| ( |2002| l, except for few differences for the values of 
the pericenter radius. However, the adopted values in our simu- 
lations make them resembling the Barnes' videos. In some cases 
we had to invert the spin of interacting galaxies to match the 
observed velocity gradient. The number of particles in each sim- 
ulation is 95040 slightly larger than in |Barnes ( 2002[ l. 

We then developed the interactive system based on Barnes' 
command "snapview" which allows us to display/rotate simu- 
lation in 3D space. We thus improved it by generating the pro- 
jected image (morphology) and velocities from which we can 
mimic the IFU observation. After matching the morphology, the 
velocity of each particle is projected to line-of-sight direction. 
Then we mimicked the observation of an IFU by calculating the 
mean velocity and velocity dispersion for those particles that are 
projected into an IFU pixel. We have first tried to reproduce the 
gross morphological features and then we have tried to catch 
the kinematics by dithering the IFU grid within half IFU pixel, 
and by rotating the system within ~ 5-10 degrees. Figure 2 dis- 
plays the final result after rotation for morphology, velocity field 
and dispersion maps which compare the best to the observations 
within the adopted grid of models. 

To measure the quality of the model, three sets of parameters 
are considered: the ones from morphology, from velocity field 
and from dispersion maps, respectively. Each set of parameters 
is graded from (failure) to 2 (good fit). For each of them the 
following criteria have been considered: 

- Morphological parameters: only large scale morphological 
structures were considered, including the presence of multi- 
ple nuclei, bars, ring, arms, helix etc.; 

- Velocity field parameters: orientation of the main velocity 
gradient(s), but not their amplitude(s); 

- Dispersion map parameters: position of the minima and max- 
ima, but not their amplitudes. 

The quaUty of a model is given by the final grade which is 
the sum of the grades for each three sets of parameters. Three 
of us have independently classified each selected model (LA, 
FH and MP). Individual grades has been then compared after 
the completion of the whole classification. An excellent agree- 
ment have been found between the individual classifications, 
with only three major disagreements within 26 objects. The final 
grade is the median value and has been provided after a face- 
to-face meeting between the three classifiers. Table 2 provides 
the individual grades as well as the resulting final one. We as- 
sume in the following that a robust model has been obtained 
when the final grade is equal or above 4, while median grade 
of each three sets of parameters has to be above or equal to 1. 
This means that we are able to reproduce simultaneously mor- 
phology and kinematics in a quite robust way. We then find that 
17 starbursts among 26 are robustly reproduced by our simple 
modelling method (see Table 2 and Figure 2). 

Let us now examine whether or not these robust fits gives 
convincing credit to the merging hypothesis. The main limita- 
tion of this exercise is obviously given by the discrete number of 
orbits and of mass ratio between the modelled encounters. Since 
we allow to invert the spin of one or two of the interacting galax- 
ies, the total number of models is 48. Notice that we had to invert 
the spin for 5 cases among the 17 robustly modeled starbursts. In 
fact our initial motivation to choose such a methodology, besides 
the obvious question of time consuming, was coming from our 
initial experience. For each individual detailed modelling (see 



Table 2. Classification of the quality of the model in reproducing 
morphology and kinematics. 



lAU name 


Class" 


CI* 




CI*™ 


CI" 




J033210.25-274819.5 


NR 


212 


5 


222 


6 


222 


6 


6 


J033212.39-274353.6 


RD 


- 


- 


_ 


_ 


_ 


_ 


6 


J033213.06-274204.8 


SR 


221 


5 


122 


5 


211 


4 


5 


J033214.97-275005.5 


NR 


101 


2 


111 


3 


111 


3 


3 


J033217.62-274257.4 


NR 


022 


4 


212 


5 


122 


5 


5 


J033219.32-274514.0 


NR 


112 


4 


202 


4 


112 


4 


4 


J033219.61-27483L0 


NR 


121 


4 


121 


4 


121 


4 


4 


J033219.68-275023.6 


RD 


- 


- 


_ 


_ 


_ 


_ 


6 


J033220.48-275 143.9 


NR 


Ill 


3 


110 


2 


010 


1 


2 


J033224.60-274428.1 


NR 


222 


6 


122 


5 


222 


6 


6 


J033225.26-274524.0 


NR 


120 


3 


111 


3 


110 


2 


3 


J033226.23-274222.8 


SR 


- 


- 


_ 


_ 


_ 


_ 


6 


J033227.07-274404.7 


NR 


222 


6 


122 


5 


222 


6 


6 


J033228.48-274826.6 


NR 


111 


3 


120 


3 


220 


4 


3 


J033230.43-275304.0 


NR 


122 


5 


211 


5 


212 


5 


5 


J033230.57-274518.2 


NR 


010 


1 


101 


2 


100 


1 


1 


J033230.78-275455.0 


RD 


- 


- 


_ 


_ 


_ 


_ 


6 


J03323L58-27412L6 


RD 


- 


- 


_ 


_ 


_ 


_ 


6 


J033232.96-274106.8 


NR 


121 


4 


122 


5 


221 


5 


5 


J033233.90-274237.9 


NR 


212 


5 


221 


5 


222 


6 


5 


J033234.04-275009.7 


SR 


211 


4 


121 


4 


121 


4 


4 


J033234. 12-273953.5 


NR 


010 


1 


120 


3 


020 


2 


2 


J033237.54-274838.9 


RD 


- 


- 


_ 


_ 


_ 


_ 


6 


J033238.60-27463L4 


RD 


- 


- 










6 


J033239.04-274132.4 


NR 


121 


4 


120 


3 


120 


3 


3 


J033239.72-275154.7 


NR 


220 


4 


220 


4 


220 


4 


4 


J033240.04-274418.6 


NR 


220 


4 


220 


4 


220 


4 


4 


J03324L88-274853.9 


SR 


222 


6 


222 


6 


222 


6 


6 


J033244.20-274733.5 


NR 


210 


3 


100 


1 


100 


1 


1 


J033245. 11-274724.0 


SR 


222 


6 


222 


6 


222 


6 


6 


J033248.28-275028.9 


SR 


221 


5 


221 


5 


221 


5 


5 


J033249.53-274630.0 


NR 


201 


3 


100 


1 


001 


1 


1 


J033250.53-274800.7 


NR 


222 


6 


122 


5 


122 


5 


5 



been improved by the intensive utilization of GPU (execution 10 times 
faster than CPU). 



Notes: 

" Morpho-kinematical classification (RD: rotating spiral disks, NR: 
non relaxed systems and SR: semi-relaxed systems, see section 1). 
'' Quality grade of the model for the three sets of parameters (morphol- 
ogy, velocity field and the cr-map). Each one is graded from 2 (good 
fit) to (failure). The last grade is the sum of the three previous ones. 
^ The final classification grade is the median value of the three grades 
attributed by each individual and ranges from 6 (excellent agreement) 
to 4 (robust fit) and then to 3 and below (not reliable fit). 



section 2.1) our initial guess for the orbit and mass ratio had been 
close to the final result, or in other words, only few changes of 
the latter quantities have allowed to recover the amplitudes of 
the kinematical parameters as well as most of the morphologi- 
cal details. It is indeed possible that the choice of the orbits by 
Barnes ( 2002[ ) is quite representative of a cosmological distribu- 
tion of merger orbits. Besides this, the use of 1:1 and 3:1 mass 
ratio might be also sufficient to reproduce the gross features of 
many major mergers. 

Our models, including the robust ones, are not unique and 
this probably applies also to detailed models that have been dis- 
cussed in section 2.1. The important question to address here 
is the possible degeneracy of the methodology applied here, 
that would be the case if the number of model parameters is 
significantly larger than the number of constraints. To test this 
we propose to consider two robust models, one of an on-going 
merger with two identified nuclei and one of a merger remnant. 
J033210. 25-274819. 5 is a galaxy showing two components with 
a strong colour difference (see Fig. 2) of 1.5 magnitude in ob- 
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served {b - z)ab (see |Neichel et al.||2008] l. To test the merger 
hypothesis we simply assumed that the two components are the 
relics of a 3:1 mass ratio merger with a DIRect orbit just before 
the first passage. This galaxy shows two peculiarities in its kine- 
matics. First the dynamical axis is seriously offset from the main 
optical axis and it points towards the secondary blue component 
and second, the dispersion peak is also offset in the direction of 
the secondary component. The constraints that are needed to be 
reproduced are: 

- the location of the secondary component (2) and the presence 
of a small bar (1) in the center of the main component; 

- the main dynamical axis (2) that is shifted towards the sec- 
ondary component as well as the presence of a secondary dy- 
namical axis (1) that follows the main optical axis of the main 
component; 

- the location of the dispersion peak (2) as well as the minima (1 
to 3) in the dispersion maps. 

The parameters required to reproduce these > 9 constraints in- 
clude only the mass ratio, the elapsed time during the merger 
and the orbit, i.e. a much smaller number of parameters when 
compared to the number of constraints. This small number of 
parameters is due to the fact that many parameters are fixed in- 
cluding the parabolic orbit, gas fraction and pericenter radius (as 



it is the case for 3 : 1 merger in Barnes ( 2002 1). Other parameters 
such as the total mass, the baryonic fraction, the profile of dark 
and baryonic matter have been also fixed during the simulation, 
consistently with the fact that we are not reproducing amplitudes 
of the kinematical properties. The second galaxy is J033232.96- 
274106.8 which is a compact galaxy. Fig. 2 shows not only the 
morphology and the kinematics but also the residual image after 
having removed the best-fitted luminosity profile that is a n=l 
Sersic index with a 6.4 HST/ACS-pixel disk radius. The residual 
shows a so-called helix structure that is reproduced by the sim- 
ulation (gas component) which also reproduces the dynamical 
axis, the structure of the velocity field as well as the dispersion 
peak and most of the minima. As for the former example, in 
this (a 3: 1 POLar) merger the number of free parameters is very 
small (3) and far below the number of constraints (8) to repro- 
duce. 

Another degeneracy that might affect the modeling process 
is the one associated with the uniqueness of the best model itself. 
Indeed, even if the number of constrains exceed the number of 
model parameters as discussed above, there remains the possi- 
bility that at constant number of free parameters, several differ- 
ent models could match the observations with a similar quaUty. 
This model-degeneracy was extensively discussed by [Barnes &| 
[Hibbard (2009). They pointed out that such a degeneracy can 
be broken by identifying specific features in the phase space (ie, 
the 6D space of positions and velocities) after the encounter, be- 
cause such features (e.g., tidal tails) allow us to trace back the 
initial configurations of the progenitors. They showed that such 
a methodology allowed them to robustly constrain "the disk ori- 
entations, viewing angles, time since pericenter, pericentric sep- 
aration, and scale factors", while they did not examine "errors in 
center-of-mass position and velocity". In spirit, our approach is 
very similar: we examined the 6D phase space using projections 
of different moments of the phase function (morphology, veloc- 
ity field, and velocity dispersion maps) and tried to reproduce 
specific signatures, which are listed above. Hence, we are quite 
confident that the methodology used in this paper allows us to 
constrain efficiently the merging phases as well as the disk in- 
clinations provided that the time since pericenter does not corre- 
spond to the latest merging phases where such specific features 
tend to vanish. But strictly speaking, it is clear that it remains 



difficult, and probably even impossible, to claim that these mod- 
els are truly unique, something which is anyway inherent to any 
modeling work, whatever the adopted methodology and/or data 
quality are. 

In the sample of 33 emission line galaxies there are six ro- 
tating spiral galaxies, five of them being isolated and one be- 
ing in interaction with an elliptical galaxy (see Fig. 1). Another 
galaxy is almost similar to a rotating disk although it experi- 
ences a satellite infall. Among the other 26 galaxies, 17 of them 
have their morphologies and kinematics robustly reproduced by 
a merger model with a number of parameters that is far below 
the number of constraints provided by the observations. There 
are nine other galaxies for which our modelling as mergers ap- 
pears less secure. All these galaxies have anomalous velocity 
fields and peculiar morphologies and generally their dynamical 
axes show significant offset with to the main optical axes. These 
galaxies have thus similar properties than those of the robustly 
modelled mergers or remnants discussed above. Some of them 
are obviously in strong interaction (e.g. J033220.48-275 143.9) 
or are very likely merger remnants (e.g. J033214. 97-275005. 5 
and J033230.57-274518.2) from their extremely distorted mor- 
phologies and kinematics. The larger uncertainty in modelling 
them could be due to the limitations of the templates used here. 
We classify them as possible mergers or merger remnants in the 
following. 

2.3. Distribution of mass ratio, merger temporal piiases and 
orbits 

In summary, we find that among 33 distant galaxies, 17 are ro- 
bustly and 9 are possibly reproduced by models of major merg- 
ers. In our simplified model, our goal is just to identify which 
configuration (phase, orbit, mass ratio, & pericenter, see Table 
3) is able to reproduce both morphology and kinematics. There 
are several possible biases in such an exercise and several of 
them have been discussed above. It is however interesting to ex- 
amine the overall distribution of the configurations parameters 
that could reproduce the distant starbursts as mergers. 

First, we may consider the mass ratio between the two inter- 
lopers. For each configuration in which the two interlopers can 
be identified, we have used z-band photometry to calculate the 
mass ratio (see Table 2). However many starbursts have been 
identified with merger remnants for which we derive the mass 
ratio from the modelling. Figure 3 (top) shows a distribution 
with two peaks at 1:1 and 3:1 (log(M2/Mi)=0 and -0.48, respec- 
tively), which are obviously related to the adopted methodol- 
ogy. Photometric estimates of the mass ratio can be done for 
mergers before the second passage (see Figures 1 and 2) when 
the two components can be separated. Nearly half of the sam- 
ple possess photometric estimates of the mass ratio, and they 
draw a smoother distribution, ranging mostly from 0.25 to 0.65 
for M2/M1 (Figure 3, middle). One may wonder how can be 
derived in such a way the mass properties of a disturbed dark 
matter component, especially for the minor interloper that is 
likely harassed during the event. Stewart | ( |2009a| l studied such 
configurations (see their Figure 2, right panel for gas-rich z=l 
galaxies), and found that the M2/M1 stellar value ranges from 
1/3 to 3/2 times the values for the dark matter, assuming stel- 
lar masses in the range of 10'" to 10" Mq, respectively. Figure 
3 (bottom) shows the distribution of dark matter ratio after ap- 
plying the correction suggested by Stewart (2009a). The main 
difference between the top and middle/bottom panels of Figure 
3 is the vanishing of the 1 : 1 peak: it is not surprising that equal 
mass mergers are rarer than 2:1 or 3:1 mergers, and indeed one 
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Fig. 2. Comparison between observations and models for the 26 distant starbursts that show non or semi-relaxed properties from 
their morphologies and kinematics. The top sketchy boxes indicates the signification of each panels. See Table 2 for the classification 
of the modelling and Table 3 for the general properties of these galaxies. 



F. Hammer et al.: The Hubble sequence: just a vestige of merger events? 



J033230.57-274518.2 (POL 1:1 ?-peri=0.4) 



J033232.96-274106.8 (POL 3:1 Vri=0-2) 





- i -il 







































3 
1 



! I 



J033233.90-274237.9 (RET 1:1 rperi=0.4) 



J033234.04-275009.7 (RET 1:1 rpen=0.2) 





J033234. 12-273953.5 (INC 1:1 rpe,,=0.4) 




J033239.04-274 132.4 (POL 1:1 ^,.,=0.4) 




4 1^ i 

Ini 



1^ 



\: * I. 



J033239.72-275 154.7 (DIR 3: 1 rpe„=0.2) Spin inv. for the small galaxy J033240.04-274418.6 (RET 3:1 Vn=0-2) 




I J. ■ ■ 




L V I. 

* L" * I. 



J033241.88-274853.9 (INC 1:1 rpe,,=0.4) 



J033244.20-274733.5 (INC 1:1 Vn=0.2) 




'1 [^■'S |] 



J033245. 11-274724.0 (INC 3:1 rpe,-,=0.4) 







' 6*1 i'^" 


E 










J033249. 53-274630.0 (DIR 1:1 rpe,,=0.2) 5pm inv. for the main galaxy J033250.53-274800.7 (POL 3:1 ^,1=0.2) 





Fig. 2. continued 



8 



F. Hammer et al.: The Hubble sequence: just a vestige of merger events? 



15 



n 1 1 \ \ \ ^ \ \ \ 1 

" * 6 rotating spirals (RD), inciuding; 

5 isoiated spirals and 1 in major interaction 
" ' 1 spiral with a satellite infall (M2/M1=1/18) 
~ * 26 NR or SR, including: 
10 - 17 robustly reproduced 
by major mergers 



I I I" 

_ 13 rrlergers with photometric! estimates of IV12/M1 



. afterLorrection from Stewart (2009) 




iJZLi L 



-2-10 1 
log(M2/Ml) 

Fig. 3. Top: distribution of the logarithm of the mass ratio, Ml 
representing the mass of the galaxy observed by IMAGES. The 
full black histogram corresponds to the 24 galaxies the nature of 
which we have robustly identified, including the 17 robust merg- 
ers, the 5 isolated rotating disks (arbitrarily set at M2/M 1=0.01), 
the spiral galaxy in interaction and the galaxy with a satellite in- 
fall. The vertical dotted lines represent the limit of major merg- 
ers (mass ratio between 5 and 1/5). Middle: same as above but 
for galaxies with separated components (see Fig. 2) for which 
we have been able to estimate the mass ratio using the z-band 
photometry from ACS, which corresponds to the rest-frame V- 
band. Bottom: same as above but after applying the correction 
by(Stewart (2009a) to recover the dark matter M2/M1 ratio. We 
have generalised this correction factor to values of M2/M1 (stel- 
lar) different than 0.3. 



can notice that a large fraction of 1 : 1 mergers are not robustly 
modelled. Both distributions are overwhelmingly dominated by 
major mergers (all but the satellite infall Puech et al. 2007 a[ l. 
The overall distribution shows the scarceness of events involv- 
ing a galaxy more massive than the observed one, since those 
are rarer due to the exponential drop of the mass function to- 
wards the massive end. The quasi absence of minor merger may 
have a different meaning because minor encounters should be 
numerous at Zmedian=0.65 (e.g. Davies et al.|2009[ l. In fact minor 
mergers are expected to affect less and in a more sporadic way, 
kinematics, morphology and star formation (see also Ho pkins et| 
|al.|2 008 and discussion in section 4.1). Overall the distribution 
of mass ratio seems consistent with a modelling of most distant 
starbursts as major mergers as shown in section 2.2. 

Figure 4 shows how the modelled galaxies are distributed 
during the various temporal phases of the merger The com- 
bination of constraints from large-scale kinematics and from 
detailed morphology generally leaves few doubts about the 
merger-phase. For example for J033224. 60-274428.1 (see Fig. 
2) the collision could not be reproduced by a second pas- 
sage, that would not fit both the morphology and the disper- 
sion peak location in direct or inclined orbits. Furthermore, 
we believe that most galaxies have their phases quite ro- 



bustly identified. This is even true for several of the nine 
"possible" mergers, including J033220.48-275143.9 for which 
the disturbed morphologies of both components evidences a 
phase between the first and the second passage. Similarly, 
if J033214.97-275005.5, J033225.26-274524.0, J033228.48- 
274826.6, J033234. 12-273953.5, and J033244.20-274733.5 are 
really mergers their highly distorted morphologies or their com- 
pactness is difficult to understand if they were not near the fusion 
stage. For J033238. 60-274631.4 (see Fig. 1, right) it is very plau- 
sible that the interaction is before a first passage, although it is 
unclear whether it could be a simple fly-by or a first stage of a 
merger. 

The result shows a relatively equal distribution of the 
merger-phases in the IMAGES sample of distant-starburst galax- 
ies. In Fig. 4, we have added the 5 rotating spirals with warm 
gaseous-disks from their low values of V/o" (see |Puech et al. 
2007a). These galaxies could well correspond to a very last 
phase: a relaxation after their disks has been rebuilt. This would 
also explain why these galaxies are forming stars efficiently, 
mostly in their outskirts (see Neichel et al.|2008) : they are still 
fed by the late infall of the gas particles that have been expelled 
at larger radii by the collision. 

Figure 4 draws an evolutionary sequence in which all dis- 
tant starbursts can be identified to a major merger phase and 
are subsequently modelled. This sequence is complementary to 
that drawn by Hammer et al. (2005 i (see their Figure 6). Notice 
that distant starbursts represent a significant fraction of distant 
galaxies as they correspond to 60% of the galaxy population at 
Zmedian=0.65 ( jHammer et al.|1997] l. 

Conversely to the mass ratio and the merger phase, an ac- 
curate determination of the orbit is much more difficult, which 
is possibly due to the adopted methodology. Indeed we may 
have lost some configurations especially for phases during or 
after the fusion. Another difficulty is a possible degeneracy be- 
tween different orbits. For example, the galaxy morphology of 
J033234. 12-273953.5 may be well reproduced by a retrograde 
merger (even better than with the adopted inclined orbit) al- 
though we have not been able to reproduce the location of the 
dispersion peak. Fig. 5 shows the distribution of orbits within 
the sample of 26 mergers or possible mergers plus the rotat- 
ing spiral in interaction for which the orbit is not constrained 
at all (labelled N/A). Only 12 galaxies have their orbits robustly 
determined, i.e. galaxy morphologies and kinematics cannot be 
reproduced by other orbits. It also shows a lack of direct and ret- 
rograde orbits, which could be also related to the methodology. 
Indeed, we have tried to reproduce both kinematics and mor- 
phological details such as starbursting bars, rings and helicoidal 
structures that are generally associated to inclined and polar or- 
bits. 

Only detailed modelling may solve the degeneracy in the or- 
bit determination and help to verify whether the orbital distribu- 
tion is consistent with the hierarchical model of galaxy forma- 
tion. In this section we have just demonstrated that for the many 
distant starbursts having complex morphologies and kinematics, 
most of them can be reproduced by a simple modelling of ma- 
jor mergers. The next step is to estimate whether this is just a 
coincidence, or if mergers may explain the elaboration of spi- 
rals, as suggested in Fig. 4. Robertson et al. (20061 and many 



subsequent studies (Ho pkins et al.|2009at|Governato et al.|2009| l 
have shown that within the conditions of sufficient gas richness 
-generally assumed to be over 40-50%- major mergers lead to 
the rebuilding of a significant disk. 
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Fig. 4. Distribution of the merger-phases for the modelled distant 
starbursts. Full black histogram corresponds to the same galax- 
ies than in Fig. 3 (robust mergers/interactions and isolated rotat- 
ing spirals). Phase 1 corresponds to the approaching phase until 
the first passage, phase 2 corresponds to the time elapsed be- 
tween the first and the second passage, and phase 3 to the time 
after the second passage and before the elaboration of the rebuilt 
disk. During phase 3 the galaxy may take the appearance of a 
chaotic morphology or a central starburst , always accompanied 
with chaotic velocity field and dispersion peaks clearly offset 
from the mass center. Phase 4 correspond to the disk rebuilding 
phase, for which the disk is detected, the rotation is evidenced 
but could be offset from the main optical axis, and the disper- 
sion peak(s) is (are) closer to the mass centre. The last phase 
corresponds to rotating spirals with regular velocity fields and 
a dispersion peak centred to the mass centre (e.g Flores et al. 
|2006l l. 



3. Gaseous content of distant starbursts 

3.1. Gas and stellar masses 

The gas content of distant galaxies is still poorly known. Large 
interferometric sub-mm baseline (ALMA) will be essential espe- 
cially in providing estimates of the molecular gas that is related 
quite intimately with star formation ( Gao & Solomon 2004 1. We 
may however make use of galaxies in the local universe, for 
which the surface densities of star formation and gas are ob- 
served to follow a Schmidt-Kennicutt la w, '^sfr ~ '^las ^ over 
more than 6 orders of magnitude in SFR (Kennicutt|1998| ). This 
empirical relation is usually explained by a model in which the 
SFR scales with density-dependent gravitational instabilities in 
the gas. 

It is plausible that Zmedian-0-65 starbursts, including those in- 
volved in mergers are following the same relation as local star- 
bursts and mergers. Gaseous masses have been estimated by 
|Puech et aL] (|2009b ) by assuming that distant starbursts followed 
the Schmidt-Kennicutt relation. It generates gas surface densi- 
ties that range from 10 to 250 Mq pc^^, i.e. intermediate be- 
tween normal spirals and ULIRGs. Such a method has been al- 



Fig. 5. Distribution of the interaction orbits adopted to model 
the 26 mergers plus one interaction in the sample of 33 distant 
starbursts observed with IMAGES. Full black histogram corre- 
sponds to the 12 objects for which we have robustly identified 
the orbit. The x-axis represents the different orbits displayed in 
the Figure. 



ready used to derive the gaseous mass fraction of very distant 
galaxies (lErb et al. 2006). Because the local estimate of the gas 
surface density is independent of the IMF, |Erb et al.| ( 2006| no 
ticed that by using the same Salpeter IMF adopted by ~~ 
( |1998| l, the derived gaseous mass is in principle an 



Kennicutt 



MF inde- 
pendent quantity. By accounting for the uncertainties of the SFR 
derived from Spitzer/MIPS 24/im measurements (Salpeter IMF 
Kennicutt 1998), Puech et al. (2009b) derived error bars ranging 
from 0.04 to 0.4 dex for gaseous mass estimates. 

Ironically, the wealth of UV to mid-IR data in IMAGES 
does not provide yet a better accuracy in estimating the stel- 
lar masses. Indeed at all wavelengths most of the emission are 
due to massive stars and not to the main sequence stars that 
constitute the bulk of the stellar mass. The results are heav- 
ily dependent on the IMF, the assumed history of star forma- 



tion and on various recipes for the extinction law. Puech et al. 



(2008) (see their Appendix A) have analysed the systematic ef- 
fects in adopting different schemes for the stellar mass estimates. 
They have adopted stellar masses {M steiiar.Bm) calculated from 
rest-frame K-band magnitudes, using the methodology of [BeUj 
|et al.| ( |2003| ) to correct for massive stars and assuming a "diet" 
Salpeter IMF. By comparing the evolution of Msteiiar/LK to that 
from comparable studies, they found that at Zmedian=0-65, the 
Bell et aL] ( 2003 i method overestimates the stellar mass by ~ 
O.I dex when compared to that derived from|BruzuaI & Chariot 



(2003)(BC03). Furthermore Mar aston et al.K2006| ) showed that 
by including TP-AGB stars, the Msieiiar/LK is overestimated by 
an additional ~ 0.14 dex. Thus stellar-mass estimates from lBelll 
et al. (2003 1 (Msteiiar.BO^) appear to be maximal estimates for the 



stellar masses of the IMAGES galaxies considered here. It also 
gives an approximation of what could be estimated from a com- 
bination of the Salpeter IMF and |Bruzual & Charlot] ( |2003] ) mod- 
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Fig. 6. Gas mass derived fromlPuech et al 



tion of the stellar mass obtained from the Bell et al. 



(|2009b| as a f unc- 
([2003]) pre- 



scriptions. The arrows in the upper left comer indicate how the 
stellar mass may be affected by different prescriptions of the 
galaxy synthesis population models (BC03 and M06 indicate the 
change by adopting Bruzual and Chariot, 2003 and Maraston, 
2006 models, respectively), or of the IMF (S55 and K02 show 
the change due to the adoption of an IMF from Salpeter, 1955 or 
from Kroupa, 2002, respectively). The full and dashed lines indi- 
cate a gas fraction of 50% and a gas fraction from local galaxies 
(HI gas, 'Schiminovich"2008), respectively. The various symbols 
correspond to the morpho-kinematic classes defined in section 1: 
rotating spirals (blue circles), heavily non-relaxed systems (red 
triangles) and semi-relaxed systems (green squares). Only error 
bars lower than 0.3 dex have been represented in this plot, and 
values with larger error bars are represented as open symbols. 



els (see arrows in Fig. 6): the IMF effect (-1-0.15 dex) is compen- 
sated by the overestim ate due to the,Bell et al. (|2003) methodol- 



ogy when compared to [Bruzual & Chariot (2003 1. Note that the 
Salpeter IMF is the maximal IMF and also violates the maxi- 
mum disk constraints for local spirals (see e.g. Bell et al., 2003). 
Adoption of other IMFs (for example, Kroupa|[2002[ ) unavoid- 
ably confirms that MsieUar,BQ'i overestimates the stellar mass (e.g. 
Bell et al. 2003 1. Note that the above applies for the whole popu- 



lation of galaxies considered here, and that the estimate of stellar 
masses for individual galaxies is still challenging. 

Figure 6 shows how the gaseous masses are distributed 
against the maximal estimates of the stellar mass (see also Table 
3). All distant starbursts but one show gas fractions intermediate 
between the local values and 50%, and there is a correlation be- 
tween their stellar and gaseous masses. Such a correlation may 
be expected as we have selected starburst galaxies on the basis of 
their [Oil] equivalent widths and a proxy of their stellar masses 
((Wo(0//) > 15A & MjiAB) < -20.3, see IMAGES-I). 



3.2. Gas fraction in distant starbursts and tlieir progenitors 

The median gas fraction of the sample is 31%+1%, which 
is much larger than the corresponding fraction for the Milky 
Way (12%, Flynn et al. 2006) or for M31 (5%, |Carignan et al. 
|2006[). It is approxi mately twice the value found by ALFALFA 
( Schiminovich|20()8 ). Different prescriptions for synthesis mod- 
els or IMF would lead to larger gas fractions in the distant sam- 
ple: for example, combining a Maraston model with a Kroupa 
IMF would divide the stellar mass by a factor ~ 2, and thus 
the median gas fraction would become 47%. It is important 
to stress that we generate minimal gas fractions by applying 
maximal stellar masses. An independent confirmation of higher 
gas fractions in distant galaxies is provided by Rodrigue s et al.| 
(j2008 ) who studied the O/H abundances in the gaseous phases 
of similar distant starbursts. They found that gas fraction should 
reach ~30% at z~ 0.65 if one assumes ~10% the gas fraction 
in present-day galaxies. Given the large uncertainties in our es- 
timates of gas fraction, an independent confirmation is certainly 
reassuring. 

Figure 7 shows how the gas fractions are distributed in the 
different merger phases defined in section 3. There is no specific 
trend between the two quantities. In the frame of a merger sce- 
nario, this is not unexpected because it results from the balance 
between two effects. Let us consider two sets of z=0.6 galaxies, 
one (hereafter called A set) including galaxies in later merger 
stage (after the fusion) and the second (hereafter called B set) 
with galaxies before the first pass. By construction, progenitors 
of A galaxies are lying at higher redshift and should have lower 
luminosities (and masses) than B galaxies, and then should have, 
on average, a significantly higher gas fraction. This is somewhat 
compensated by the fact that A galaxies have had a supplemen- 
tary star formation that is induced during the merger, which has 
transformed a part of the gas into stars, reducing the gas frac- 
tion. Given the small statistics (5 to 7 galaxies in each phase), 
we believe that Figure 7 is consistent with the merger hypothe- 
sis. Indeed the main question is whether the progenitors of the 
observed galaxies are sufficiently gas-rich to lead to disk rebuild- 
ing. 

There are two possible means to estimate the gas fraction in 
merger progenitors. Hammer et al. ( 2009b) have modelled the 
stellar populations from deep spectroscopy of a merger remnant 
assumed to be in a disk rebuilding phase. They found that ~ 50% 
of the observed stellar mass had been formed during the merger 
by comparing the stellar population ages to the merger dynam- 
ical time. This implies that the progenitors of this system were 
on average more than 50% gas-rich, supporting the evidence that 
this system is rebuilding its disk after a major merger Here we 
use the characteristic doubling time, r5/rR=Mj,f//a,-^03/SFR, to 
estimate the stellar mass that has been formed during the event, 
and thus the gas fraction in their assumed progenitors. Assuming 
that our maximal stellar-masses are a good representation of the 
stellar mass for a Salpeter-IMF and a Bruzual & Chariot (2003 1 
synthesis model, this quantity is IMF independent, as our SFR 
estimates have been done using the Kennicutt ( 1998 ) calibration 
(Salpeter-IMF) to the IR luminosities ( Puech et al.||2009b) see 
details). 

Figure 8 gives the distribution of the characteristic stellar 
mass doubling times for the 34 IMAGES starbursts. It is re- 
markable that their median value takes its minimum near the 
fusion, which is expected in every model of mergers. We have 
used the models shown in section 3 to estimate the time each 
galaxy spends in each of the phases (see time values in Table 
3). It assumes a rotating period time of fr„,= 1.2 x 0.25 Gyr for 
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merger phase merger phase 




Fig. 7. Gas fraction distribution as a function of the merger 
phases. The former corresponds to a minimal value of the gas 
fraction in the sample of IMAGES starbursts (see text). The var- 
ious symbols correspond to the morpho-kinematic classes, as in 
Fig. 6. Note the presence of the massive dry merger at phase=2 
(see the inserted image in Fig. 6) with a very low gas fraction 
( Yang et al.| [2009). The full blue line shows the median gas- 
fraction in the progenitors derived from the gas consumption 
during the merger (see text). The dot-long-dashed line shows the 
same, but for a combination of [Maraston et al.]p006| l models 
and Kroupa IMF. The black solid line gives the median of the 
gas fraction values. 



a galaxy with the mass of the Milky Way (see Barnes 2002 1. 
We then scale the merger time with the observed baryonic mass 
(assuming MMw.barvonic-5.5 10'" Mq) and a lso apply a correc - 
tion for the merger mass ratio as described in 'Jiang et al. ( 2008| l. 
Besides this, we calculate for each phase the median Ts fr that 
is considered to be its effective star formation time. For a given 
starburst assumed to be in a given phase, we may calculate the 
fraction of gas that has been transformed into stars during each 
previous phases of the merger. Table 3 (column 10) gives the 
time spent by each galaxy in the previous merger phases, that, 
after combination with the median Ts fr, provides us with an es- 
timate of the gas mass that has been transformed into stars during 
the merger. 

Figure 7 also shows the resulting (median) distribution of 
the progenitor gas-fractions (see also Table 3). In 75% of the 
progenitors, the gas-fraction is in excess of 40% or of 58% if we 
adopt a Maraston-Kroupa combination. The high gas fraction in 
progenitors is robust because: 

- it only requires a modest star formation in the merger to reach 
large values since we already find a high fraction of gas (median 
31%) in the observed starbursts ; 

- the gas-to-star transformed mass during the merger is inde- 
pendent of the observed mass and thus of the IMF choice, since 
both Ts FR and the merger time scale with mass. 



Fig. 8. Characteristic stellar mass doubling times as a function of 
the merger phases, as they are described in Fig. 4. The various 
symbols corresponds to morpho-kinematics classes as in Fig. 6. 
The full blue line indicates their median value and the dashed 
red line indicate the time during which the observed starbursts 
are involved in the merger (median value for each phases 1 to 4), 
as obtained from the simulations. It slightly decreases at phase 4 
because the masses of the 4 rebuilding disks (median value, 2.1 
10'" Mq) are smaller than the mass of the 9 galaxies in phase 3 
(median value, 6.3 10'" Mq). 



4. Discussion 

4.1. A self-consistent explanation of the galaxy 
transformation during the last 6 Gyrs 

In section 2 we argued that distant starbursts have morphologies 
and kinematics consistent with major mergers or their remnants. 
In section 3 we show that they have large gas fractions, and that 
their progenitors would have to have gas fractions above 40- 
50% to account for the stellar mass produced during the merger. 
Therefore most starbursts at Zmedian=0-^5 - those with anoma- 
lous morphologies and kinematics - are consistent with gas-rich 
merger phases leading to rebuilt disks. 

Our interpretation of the morpho-kinematic evolution (see 
Table 1) is then straightforward: ~ 6 Gyrs ago, 46% of the 
galaxy population was involved in major mergers and most of 
them (75% x 46%= 35%) were sufficiently gas rich to rebuild 
a disk. Those can be considered as progenitors of the present- 
day numerous spirals - although this deserves a careful analysis 
of the exchanges of angular momenta - while the others could 
be progenitors of E/SO and of the scarce population of massive 
irregulars at the present-epoch (~ 10%, see Delgado- S errano | 
et al. 2009). Thus as much as half of the present-day spirals 
are coming from disk rebuilding from recent mergers, the other 
half being already assembled into quiescent or warm disks at 
z,nedian=0.65 (Table 1). 

More statistics are certainly needed to obtain a more precise 
estimate of the amount of gas that has been consumed during the 
different merger phases. The median time spent in each merger 
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( |Rodrigues et al.|[2008[ ). It is also consistent with the observa- 
tions of galaxy pairs at z ~ 0.6. The most robust estimate is that 



Fig. 9. b+v, i and z combined images of three additional star- 
bursts whose kinematics has not been detected due to spatial res- 
olution. The galaxy in the middle panel show a ring and possibly 
the nucleus of the secondary interloper Kinematics are needed 
to confirm it as well as to verify the nature of the two other star- 
bursts that could be in close interaction and just at the nuclei 
fusion, respectively. 



phase ranges from 0.5 to 1.4 Gyr (see Fig. 8 and also Table 
3): the scenario naturally explains why distant starbursts show 
a so important contribution of intermediate-age stars revealed 
by their very large Balmer absorption lines in their spectra (e.g. 
Hammer et al. 1997; Marcillac et al. 2006, see also Poggianti et 



al. 1999 for another perspective in galaxy-cluster environments). 



The median baryonic mass of the sample is 0.75 times that of 
the Milky Way. Their progenitors should be galaxies at larger 
redshifts, approximately 1 Gyr earlier, i.e. at z~ 0.83. At such 
redshifts the large gas fractions in progenitors is not exceptional. 
Accounting for the gas consumed during the merger, the median 
stellar mass and gas fraction of their progenitors are 7.5 10^ M© 
and 50%, respectively. In present-day galaxies within this mass 
range, the gas fraction averages to ~26% for local galaxies (from 
Schiminovich 200 8 j, and it could be not exceptional that 7 Gyrs 
ago such galaxies had twice their present gas content. 

Improvements are also required to estimate the stellar masses 
since a proxy (absolute J-band magnitude) of the stellar mass has 
been used in this study to select our sample. Combination of re- 
alistic stellar population with different ages and metal content 
has to be performed on both the whole spectral energy distri- 
bution (from UV to near-IR) and the spectroscopic absorption 
lines (Lick indices). Nevertheless, we do find that all distant 
starbursts are consistent with major merger phases, and these 
sources are strong emitters in near-IR. It is unclear whether we 
may have missed a significant population of massive starbursts 
without strong emission in near-IR. Besides this, technical lim- 
itations (see IMAGES-I) have prevented us from measuring the 
kinematics of 3 starbursts, because their emission line region are 
too compact (see Fig. 9). Their optical morphologies are also 
consistent with mergers (see Fig. 9 and its caption). 

Figure 4 - see also Figures 1 and 2 - suggests that all dis- 
tant starbursts are part of the duty cycle of the disk rebuilding 
scenario (see Hammer et al. 2005 their Figure 6). It also in- 



cludes the distant rotating spirals with high star formation and 
warm disks, which are the natural last phase of such a gas-rich 
merger event. During the elapsed time to z=0, they may trans- 
form their gas-into-stars and simply relax to form the present- 
day thin disks. During such phase, models predict that almost all 
the gas has reached the disk, which is confirmed by the fact that 
they lie on the same baryonic Tully Fisher relation as local disks 



( Puech et al. 2009b I . In passing, this brings a simple explanation 
of why some regular disks are LIRGs (Melbourne,2006) . 

The spiral disk rebuilding scenario explains the changes ob- 
served in the galaxy population, the density evolution of star 
formation and stellar mass, the dispersion of the evolved Tully 
Fisher relation ( |Flores et~aL| 2006 , Pu ech et al.|2009bl ) as well as 
the strong evolution of the metal abundances of their gas phases 



5±1% of galaxies are in pairs at that epoch (see e.g. Bell et al. 
2006, for a review): in our sample we identify only 4 galaxies 
which could have been identified as well separated pairs. It cor- 
responds to a pair fraction of 4/33x0.6= 7+3% (see Table 1 in 
which 0.6 is the fraction of starburst galaxies) in the Zmedkm=0.65 
population. Note that among these pairs, three are in phase 2, i.e. 
between the first and the second passage. Last but not least, ma- 
jor mergers may explain why the ionised gas radius is larger than 
the optical radii, especially for the starbursts with small optical 
radius (Puech et al. 2009b). Indeed, during a close encountering 
and nuclei fusion phases, the gas is heavily shocked by the colli- 
sion and could be ionised this way (Puech et al. 2009a i. Another 
alternative is that dust-enshrouded clumps may ionise the gas 
while they cannot be detected at visible wavelengths; this phe- 
nomenon has been identified in a compact dust-enshrouded disk 
(see Hammer et al.^2009b) . 

Cosmological simulations confirm the importance of the gas 
in mergers. By accounting for the gas, ( [Stewart et al. |[2009b| l 
found that almost all galaxies may have experienced a merger 
since z=2. A comparison to IMAGES observations would be 
highly desirable to sharpen this prediction. Those may be consis- 
tent with the whole formation of the Hubble sequence: a forma- 
tion of spirals by the numerous gas-rich mergers and of massive 
ellipticals by gas-poor mergers. It is now time to study the for- 
mation (or re-formation) of the present-day spiral galaxies and 
of their substructures and to catch a glimpse on how the whole 
Hubble sequence may have been formed. 

4.2. How other mechanisms contribute to the galaxy 
transformation? 

Numerous studies have attempted to describe the above evolu- 
tion by assuming different mechanisms. It has been argued that 
only a small fraction of the star formation is triggered by major 
interactions ( [Robaina et al.||2009[ [Jogee et aL||2 009). Although 
these studies are based on a large number of objects, they do not 
possess the kinematical information that is crucial to evaluate 

( |2008l l 



the presence or absence of merging. Indeed Neichel et al. 



found that most massive galaxies showing irregularities or com- 
pactness have also anomalous velocity fields. Other mechanisms 
than merging have then to explain not only morphological irreg- 
ularities but also their non-relaxed kinematics. 



We have examined the spectra (see e.g. Rodrigues et al. 
2008) of 20 IMAGES starbursts to investigate whether outflows 
can be detected. An outflow may lead to significant differences 
between the velocity of the emission-line system and that of 
the absorption-line system ( jHeckman & Lehnert |2000| . We do 
find diff'erences in only 3 (J033214.97-275005.5, J033224.60- 
274428.1 and J033225.26-274524.0) within 19 objects, at the 
level of ~ 100 km/s. Thus stellar feedback mechanisms are un- 
likely to considerably affect distant starbursts, which is under- 
standable as all of them are relatively massive with baryonic 
masses in excess of 10'° Mq. The absence of mmor mergers m 
our sample (only one case with 18:1 mass ratio) could be sim- 
ply explained by their considerably lower efficiency. Due to their 
lower impact and their longer duration ( |Jiang et ar|2 008 ), they 
are considerably less efficient to activate a starburst, and to dis- 
tort morphologies and kinematics or then they do it in a some- 
what sporadic way (see [Hopkins et al.|[2008) l. To explain both 
the stellar mass assembly and the peculiar morpho-kinematics 
would probably need an extremely large (non-observed) rate of 
minor mergers, that is certainly not consistent with the large an- 
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gular momenta of present-day disks (see e.g. Mailer & Dekel, 
2002). 



5. Conclusion: a scenario to explain the elaboration 
of the Hubble sequence? 

In this paper we considered the possibility that the formation 
of the Hubble sequence relies to a large extent on past merger 
events. We used a sample of objects around z=0.65 for which 
we have both morphology from the HST and high quality kine- 
matics (both velocity field and 2D velocity dispersion maps) 
and compared them with simulation results varying the view- 
ing angles to obtain the best fits. Although we can of course 
not prove that the origin of these objects is indeed a merger, we 
can safely say that their observed properties are well compati- 
ble with them being merger or their remnants. We want to stress 
that this result was reached using all available data, i.e. our com- 
parisons included morphology, mean velocities and dispersions. 
This strengthens our proposal very considerably. 

A merger origin of the Hubble sequence is very intimately 
linked to the disk rebuilding scenario, which has been succes- 
fully argued both from the observational point of view and from 
simulations ( Barnes|2002 Springel & Hernquist |2005 Hopkins 



|et a l. 2009a). Our quantitative estimates of gas fractions argue 
that more than a third of the galaxy population about 6 Gyrs ago 
was sufficienty gas rich to rebuild a disk after the merger Thus 
our work argues that a merger origin of the Hubble sequence, 
although it has not yet been proven, is a plausible alternative or 
channel for the formation of large disks in grand-design spirals. 
Half of the present-day spirals being in merger phases at 
naturally implies that most and probably all were 
shaped during gaseous-rich mergers at earlier epochs. These 
mergers at Zmedian-0-65 have generally begun 1 Gyr earlier, i.e. 
at z= 0.835. We may expect, that a similar amount of mergers 
had occurred a further 1 Gyr ago (from z=0.835 to z=1.07): then 
almost all spirals may have been rebuilt from their last major 
merger, during the last 8 Gyrs. This may apply to M31 ( [Hammer] 
et al.|2007 ), but the Milky Way appears quite exceptional as its 



properties imply a last merger only at much earlier epochs (10 
to 11 Gyrs ago). The Milky Way quiescent history is well illus- 
trated by its exceptional pristine halo; combined with its lack 
of angular momentum and stellar mass, this may simply indicate 
that our Galaxy has exceptionally avoided any major merger dur- 
ing a large fraction of the Hubble time ( Hammer et al. 2007) . 

A considerable task is thus awaiting us. We have to relate the 
distant starbursts to local galaxies by modelling in detail all the 
distant galaxies for which we possess detailed morphologies and 
kinematics, i.e. about 100 galaxies. Although large, such a num- 
ber is barely sufficient to describe the large variety of merger 
configurations. With such a modelling we will be able to derive 
the end-up properties of each starburst, by modelling their evolu- 
tion 6 Gyrs later, and verify whether they are consistent with the 
present-day distribution of galaxies within the Hubble sequence. 

In order to make better comparisons between observations 
and simulations it would be useful to have a larger library of sim- 
ulations covering longer evolution times and larger gas fractions. 
A better coverage of the parameter space would also be very use- 
ful, although from our experience, it will take some time before 
being able to recover detailed modelling of each of the observed 
galaxies. Such a library (e.g. 'GALMER', see |Di Matteo et aT 



|2008 ') may be useful as it could be adapted to the IMAGES ob 
servational parameters. 

Nevertheless, it was possible to reach a number of conclu- 
sions. It is instructive to note how frequent structures such as 



Fig. 10. Gas distribution in merger remnants from Barnes (2002 1 



for 3:1 mergers, and for four different orbits, at the end of the 
simulations. 



bars and rings are. We observe 9 bars and 6 rings in our sample 
of 34 galaxies, and they have colours consistent with young or 
intermediate-age stars, and as such they could have been formed 
during the merger. Figure 10 evidences that such structures are 
still persistent at late merger phases. The fraction of barred 



galaxies that we find is compatible with that found by |Sheth et al 
I ( |2008| l for the same redshift range. Bars and rings are also ob- 
served in local galaxies and may well be triggered from the ear- 
lier interactions and merging. For example, it is known that inter- 
actions can trigger bar formation (|Noguchi|1987[|Gerin, Combes] 
& Athanassoula||1990[ [Steinmetz & Navarro||1999| l. Also polar 



encounters can create rings, pseudo-rings, or spirals, which have 
characteristics similar to the observed ones (see also |Berentzen| 
et al. 2003). Finally spiral patterns can also be generated by en- 
counters, as has been shown by both simulations and observa- 
tions (e.g. |Toomre|1981||Goldreich et al.|1978|[T979HDonner,| 



'Engstrom & Sundelius 1991 Tutukov & Federova 2006 and ref- 
erences therein). Adjustment of the models is clearly a central 
issue. The helicoidal structure found in many distant starbursts 
(see examples in Figure 2) is likely due to the central torque 
described by [Hopkins et al.| ( |2009a l, and it may regulate the an- 
gular momentum transfer This structure is preponderant in the 
nuclei fusion phase and seems present in galaxies in later phases 
( [Hammer et a l. 2009b). Its efficiency in regulating the bulge- 
to-disk ratio is likely considerable (see [Hopkins et al.[[2009a| l 
although larger statistics are mandatory to verify this prediction 
and its actual role. 

Acknowledgements. We ai'e especially indebted to Josh Barnes for making pub- 
licly his simulations of galaxy mergers on his web pages as well as letting us 
using the ZENO code We are very grateful to the referee whose suggestions 
have considerably helped us improve the manuscript. 
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